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We consider a system of classical Heisenberg spins on a cubic lattice in dimensions three or more, 
interacting via the dipole-dipole interaction. We prove that at low enough temperature the system 
displays orientational long range order, as expected by spin wave theory. The proof is based on 
reflection positivity methods. In particular, we demonstrate a previously unproven conjecture on 
the dispersion relation of the spin waves, first proposed by Frohlich and Spencer, which allows 
one to apply infrared bounds for estimating the long distance behavior of the spin-spin correlation 
functions. 

I. INTRODUCTION 

Recent advances in film growth techniques, in the control of spin-spin interactions and in the 
ability to characterize magnetic materials have revived interest in the low temperature physics 
of magnetic systems. Both experimental and theoretical studies have revealed several unusual 
properties of magnetic films, such as spontaneous formation of striped patterns, reorientation 
transitions (in temperature and in the sample thickness), increase of the static magnetization 
with increasing temperature, just to mention a few [7]. It is believed that an essential role in 
determining the nature and morphology of the ordered state is played by the dipolar interac- 
tion. Unfortunately, its long-range nature and its anisotropic character make many standard 
theoretical methods and numerical algorithms inapplicable. It is therefore not surprising that, 
for instance, existence of long range order in 2D continuous spin systems interacting via a pure 
dipole-dipole interaction at low temperatures is still a matter of discussion, even at a heuristic 
level. In fact, the case of 2D lattice dipoles is a paradigmatic example of a system where the 
Mcrmin- Wagner theorem cannot be applied, spin-wave theory does not provide any resolutive 
answer and numerical simulations are difficult because of the slow relaxation dynamics associ- 
ated with the long-range nature of the interaction. New and more sophisticated methods, such 
as renormalized spin wave theory [5, 6] or block spin reflection positivity [13], are required to 
deal with this class of systems. 

At a rigorous level, even widely accepted results, such as the existence of long range orienta- 
tional order in three dimensional lattice dipole systems or the existence of the infinite volume 
Gibbs state for any given domain shape are yet to be fully proved. Many fundamental contri- 
butions to the rigorous theory of lattice dipole systems (and more generally of dipole gases on 
the lattice or in the continuum) date back to the 80 's. The use of several different techniques, 
such as reflection positivity, correlation inequalities, cluster expansion, renormalization group, 
allowed people to rigorously prove, e.g., no-screening theorems at any activity and temperature 
[10, 15], analyticity of the pressure at small activities [1, 2, 4, 12], existence of a scaling limit 
(Gaussian free field) for lattice dipoles at small activities [14] and the existence of the thermo- 
dynamic limit in the continuum in three or more dimensions [9], see [3] for a comprehensive 
review of these results. In a seminal paper, Frohlich and Spencer [10], among several other 
results, proved existence of long range order for a system of discrete lattice dipoles in two or 
more dimensions. They also described a proof for the physically relevant case of continuous 
dipoles on the cubic lattice in three or more dimensions. However, in the case of infinite-range 
interaction, they could not give a complete proof, and their argument was based on an unproven 
conjecture on the dispersion relation of the spin waves. In this paper we will give a complete 
proof of long range order for classical dipoles on a cubic lattice in three or more dimensions, 
proving in particular the aforementioned Frohlich-Spencer conjecture. The proof is based on 
reflection positivity and it extends ideas proposed in [10]. 
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II. MAIN RESULTS 

Let A be a periodic box in if , viewed as the restriction of a periodic box in M d to if . We 
assume that A is of side 2L, with L even, and we write: A = \—L + 1, . . . , L] d . We consider the 
following Hamiltonian: 

d 

Ha=J2 E SiW ij (^-y)Sl (2.1) 

*>j=l x,y£A 

where S x is a unit vector and {Sx}i=i,...,d are its components. Moreover, denoting the Yukawa 
potential by 

r Air pikx 

^(x) = / . (2-2) 



(27r) d k 2 + e 2 



the interaction matrix W(x) is defined as: 



w«(x) = ^ (-a 4 a,)r £A (x + 2ni) , x ^ o (2.3) 

and W,j(0) = ^2 a ^ (—didj)Y EA (2nL). The parameter £a is an infrared regulator that is sent 
to zero in the thermodynamic limit, i.e., lim|A|^oo £a = 0. 

For any fixed A and (3 > 0, let us denote by {■) p A the Gibbs expectation given by the 

probability measure Zp\ IIxgA d^(S x )e~ /3HA , with dfi(S) the uniform measure on the unit 

sphere and Zp t \ the obvious normalization factor. For any x <G A, given the unit vector S x , we 
define: 

< - (-ir +x *Sl (2.4) 
and denote by <r x the unit vector with components a^. Our main result is the following. 

Theorem 1 (Existence of orientational Long Range Order). If d > 3, there 
exists (3d> such that, if [3 > [3d, in the thermodynamic limit, 

IaSoo W S ^ ' &V) ^ A ~ > ' (2 ' 5) 

l_>0C 1 1 x,y£A 

with Cd{(3) a suitable positive function, vanishing at (3 — (3d- 

Using the methods of [9], it can be proved that the state {■) p A admits a thermody- 
namic limit, which we will denote by {■) p. The theorem above implies that, for (3 > {3d, the 
infinite volume Gibbs state (•) /3 is not an extremal Gibbs state. From the assumed symmetry 
of (-)p A under exchanges of the coordinate axes, and by the general theory of decomposition 
into extremal states, it follows that (-)p is a mixture of at least 2d extremal Gibbs states (pure 

phases), (-}p X \ which break rotational invariance and are characterized by 

(St)? = (-l) x+ ^l , (2-6) 

where {v\ : A = 1, . . . , 2d, . . .} are vectors obtained from some vector vo e M. d by applying 
arbitrary rotations around the origin which leave the unit cube centered at the origin invariant. 

The rest of the paper will be devoted to the proof of Theorem 1. As mentioned in the 
introduction, we will follow the same strategy proposed by Frohlich and Spencer in [10] and we 
will extend their reflection positivity method to prove, in particular, their Conjecture 7.9 [10]. 
In order to present a self-contained proof, we shall reproduce below some of the statements 
already proved in [10], including a construction of the ground states of (2.1). 
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III. PROOF OF THEOREM 1 
A. Reflection Positivity. 

Let us recall the notion of reflection positivity, adapted to the present case. If S is a unit 
vector, define 

(RiS) j = (-l) 1 ' 6 ^ S j . (3.7) 

Let TTi be a pair of planes perpendicular to the i-th direction, midway in between two lattice 
planes and bisecting A into two pieces A + and A_ of equal size. Let n denote reflection of sites 
with respect to 7r,. Clearly r^A_ = A + . We define 

{9iS) x — RiS riX (3-8) 

Let S± = {(S x }xga ± • If A is a function of S + we set 

(M)(5_) = A({(eJ) x } xeA _) (3.9) 

We shall say that the expectation (•)„ A is reflection positive (RP) iff, for an arbitrary function 
A of S+ 



A(S-)A(S + ))p A > (3.10) 



for i = 1, . . . , d. As discussed in [8] and in [10], a sufficient condition for (3.10) to hold is the 
following. Let p : M. d — > M. d be an arbitrary M. d valued function of M d . Assume that 

d 

-£ y xi>o p'(x)(M m (y)(-9^ m )y £A (x- y )>o (3.ii) 

l,m— 1 y^<o 

for all p(x) and i = 1, . . . , d. Then (-^ A is RP. 

In our context the proof of (3.11) proceeds as follows. For definiteness, let us assume that 
i = l. Let x\ > and let us rewrite 

Y £ (x) = f — e — = - f , dk± e^-^e-Myff+Z (3 12) 

Y£aW 7(2^+4 2(2^ J VkT^I 

where in the last expression = (fc 2 , . . . , fed) and similarly for xj_. If xi > yi, given any two 
functions p\ (x) , p 2 (x) : 

£ pi(x)p 2 "(y)9^ m y £A (x-y) = ^Jp- I I ^ e *x-(«x-yx) e -(* 1 -»0 V ^+5 , 



^(xj^/ki + 4 - ikxcftxjj ^(y)Vkl + 4 - ik ± ^-(y)J (3.13) 

where g± = {p{, . . . , pf ) and 02 = • • • , pa- 
using this expression we find that, if i = 1, (3.11) can be rewritten as 

1 f . f , f dk 



2(2 7 r) d - 1 y xi>0 y yi>0 y ykyn 



rf X / dy / ^ e ikx-(x ± - y± ) e -(x 1+i;i )VklT4, 



A 



■(V(x) A /ki+4-ikxe ± (x)) ^ 1 (y) v /ki+4 + zk ± ^(y)) (3.14) 

that is clearly non-negative. Let us remark that condition (3.11) is equivalent to the statement 
that Ha = H\(S-,S + ) can be rewritten in the following form [8]: 

H A (S-,S+) = H+(S+) + 9iH+(S-) - I dp(d) i C q (S_)C q (S+) , (3.15) 
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for a suitable positive measure dp. In our case, 



h+(s+)=J2 E ^(x- y )S y . 

i,j=l x,y£A + 



(3.16) 



Moreover, if i = 1 and A + = [1, . . . , L] x [— L + 1, . . . , L] d x , q is a (d — l)-dimensional vector 
and, denning = (S^, . . . , S*), 



x£A_ 



dp(q) = 



dq e-V^ 



(3.17) 



If i > 1 and/or A + is different from [1, . . . , L] x [— L + 1, . . . , L] d 1 , analogous expressions for 
C q and dp(q) will be valid. 



B. Ground states. 



Let us now show how reflection positivity allows us to construct the ground states of (2.1). 
The key remark is that the positive measure dp(q) in (3.15) induces the definition of a scalar 
product between spin configurations in A + . In particular, combining the Cauchy-Schwarz 
inequality and the inequality of arithmetic and geometric means, we find that 



/ 



d/9(q)0 i C q ($_)C q (S + ) < 



< 



dp(q) e i Cq(0iS + )C^(S + ) 



1/2 



dp(q) 0iC q (S_)C q (0i,S_) 



1/2 



< 



J d P (q) e i a l (e i s + )a i (s + ) + Jdp(q) BiC^s-)c^0iS-) . (3.18) 

If we insert this estimate in (3.15), we find that 

H A (S-,S+) > ^H A (6 t S+,S+) + ±H A (S-,0iS-) , (3-19) 

that is, either {#j<S + ,<S + } or {<S_,#j<S_} has lower energy than {<S_,<S + }. If we keep re- 
flecting in different planes, using the chessboard estimate (see Theorem 4.1 in [8]), we find 
that the energy H A (S) of a generic configuration of spins S = {>Sx}x6A is bounded be- 
low by | A| _1 X)x gA H A (S Xo ), where the spin at x in the configuration <S Xo is given by 

(— l) x (— l) Xi Sx , i — 1, ■ • ■ ,d. We now show that H A (S Xo ) is independent of S Xo . Note that 
this is not apriori obvious, since the Hamiltonian is not invariant under global rotations. Let 
S x „ = (sin # cos 0, sin6>sin</?,cos#). Then 



H A (S Xo ) 



E (-D 



xG A 
y£Z d \x 



J2(-ir- y *(Sl Q ) 2 d?Y EA ( X -y) + 



(3.20) 



• E 



-i)wsj si a i a J -r eA (x-y) 



In (3.14) the summations over x and x y factor, both for the term in the first line and for 
the one in the second line. In particular the term in the second line is zero after summation. 
So we are left with 



H A (S X0 ) = e |A| , e = £(-l) x+Xl (-8?)r eA (x) 



(3.21) 



where we used that ^(S^) 2 = 1 an d tnc rotation symmetry of F £A (x). 
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C. Infrared bounds. 



In this section we will extend the ideas used above to construct the ground states, and we 
will derive lower bounds on the Fourier transform of (tr x • a y ) „ A (also known as infrared bounds 
[11]). This will conclude the proof of Theorem 1. The main ingredients that we will need are 
Gaussian domination [11] and a refinement of the estimates on the dispersion relation of the 
spin waves, including a proof of Conjecture 7.9 in [10]. 

As a first step we map the spin system in (2.1) onto a ferromagnetic spin system via the 
mapping S x < — xr x defined by (2.4). In terms of the ct's, the Hamiltonian can be rewritten as 

d 

ff A = E E "iW&(x,y)4 (3.22) 

i,j=l x.yGA 

with W£(x,y) = (-l) x - y (-l)^ + *Wy(x-y). Note that W'(x,y) is not translation invariant. 
The discussion above shows that the ground states of H' A are spin configurations with spins all 
pointing in the same direction in space. Moreover H' A is reflection positive with respect to the 
ferromagnetic reflection (#-<?) x = oy iX . We rewrite H' A in the form 



h ' a = ~\ E E - 4)w^, y )(ai - 4) + E E «Y, w '^y) ( 3 - 23 ) 

i,j=l x.yGA i,j=l xeA y 

Note that J2 y W?-(x,y) = e l, so that 

3 

2 



H '^ = -\il E (<-4)K(^y)^i-4) + e o\M (3.24) 



*>j=l x,yeA 

Let us now define: 



*>j=l x.yGA 



J^,A(K) = <«*p{f £ E K-4- ft x + W,( x 'y)K- CT y- /i x + ^ y )})o,A (3.25) 



for /i x , x e A, real vectors. The chessboard estimate (see Theorem 4.1 in [8]) shows that 
Kp,A{E) < -ST/3,a(0) = Kfi. a (Gaussian domination). This implies ■^Kp t \{Xn)\^_ < 0, that is 

d d 

/3(IE E ( h i-K)K^y)^i-4)\ 2 Yp,A<- £ E (^-^r^^y)^-^) 

i,j=l x,y£A i,j=l x.yGA 

(3.26) 

where (-)^ A is the average with statistical weight Z^ A e~^ HA . Note that (3.26) holds apriori 
for real h x , but it extends to complex vectors [8]. In terms of the original spins, (3.26) reads 

d d 

2/3(1 E E ^-(x-y)^- eo ^4-s x |V< E E(^r^(x-y)/4-eo£i/Jxi 2 

i,j = lx,y£A x£A i,j=lx.y£A x£A 

(3.27) 

for some new vectors h x (that we are denoting by the same symbol for simplicity). Let p e 
^ L = {if m ' — m iT m 2,m 3 < 2L} and let Wjj(p) = X)xeA e lpx Wy :(x) denote the Fourier 
transform of W(x). We shall also define 

= uTi7^ E eiP ^x , Q«(P) ^ (^p^),, A • (3-28) 
' ' xeA 

Given p such that W°(p) = W(p) — e > as a matrix and choosing h x = | A| — 1 / 2 e — * px t7 in 
(3.27), we get: 

0<^ (p)0(p)^ (p)<^iy°(p), (3.29) 
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in the sense of an inequality between non-negative matrices. If, moreover, W°(p) > 0, then 

< Q( P ) < ^^(Pr 1 • (3-30) 

Eq.(3.29)-(3.30) are the key bounds. In order to make use of them, we need to study some 
properties of W°(p), in particular we need to: (i) show that W° > 0; (ii) determine the set 
S of momenta where VF°(p) has a vanishing eigenvalue; (iii) determine the behavior of its 
eigenvalues close to the zeros. Property (i) was prover in [10] (we shall reproduce the proof 
below). Moreover, in [10] it was conjectured (see Conjecture 7.9 in [10]) that, if ire is the vector 
with components (^e) 3 = 7r(l — 8t t j), then S = {ir^ : £ = 1, . . . , d} and the eigenvalue A^(p) 
vanishing at 717, satisfies A^(p) > c|p — 7T£| 2 close to tt£, for some c > 0. Our next goal will be 
to prove this conjecture. 

Let us start with showing that M /0 (p) > 0, for all p's. This is equivalent to the claim that, 
for any v £ M. d and for h x — e lpx v, 

d 

E E {KfW^-y)h{>Q. (3.31) 

i,j=l x.yGA 

By the chessboard estimate, the left hand side of (3.31) is bounded below by 
|A| _1 ^ xoeA (i?A(Wx ) — eo|w| 2 ), where the spin at x in the configuration 7i Xo is given by 
(— l) x (— l) Xi /i4 o , i = 1, ...,<£ Now note that H\(H Xa ) = e , so the proof of (3.31) is con- 
cluded. Note that if p e S then W"°(p) = 0. We now want to get a more refined bound from 
below on W°(p). In order to do this, we use the assumption that L is even (so that the side of 
A is divisible by 4) and we repeatedly reflect the left hand side of (3.31) in planes -k\ bisecting 
the horizontal bonds of the form {(2m — 1, xj_), (2m, xj_)}. Finally we repeatedly reflect in all 
possible planes 7T2, . . . , ltd- The result is a bound from below of the form: 

d 1 d 

E v l W°( P y > m E E - y)4 eoH 2 , (3.32) 



; - 1 i.j- ' \.\ \ 



with 



where: 



ul = e ^/i(^i)(_i)x+^i Vl ; 

v* t = e i ^ x ^(-l) x+Xl+x *fo(x 1 )v i , i>l, (3.33) 



fo(x) 



( 1 if x = 4fc 

1 if x = 4k + l 

-1 if x = 4k + 2 ' 

-1 if x = 4k + 3 



h(x) = f (x-l) (3.34) 



and k e Z. Given a lattice function g(x) of a single variable with period 4, we shall write 
g(x) = {g(Q),g(l),g(2),g(3)}. Note that with this convention fo(x) = {1,1,-1,-1} and 
/i(a;) = {-1,1,1,-1}. 

We now turn to the computation of the r.h.s. of (3.32). First of all, note that 
y GA( u x)* 1 ^j ( x — y) u y = if i 7^ j. This is because the double summation can be rewritten 
in the form X] x gA SyGz d \x F ij ( x ' x _ 50 > for some function F tj (x, x - y) that is odd in (x 2 - V2) 
and/or in (x 3 — y 3 ), . . . , (x d — yd), depending on the specific matrix element. Let us consider 
y eA( u x)*W /r ii( x ~ y) w y - This double summation is equal to 

u iE E e- i ^(/i(^)-/i(!/i))(_i)x-y+(^-!/i)(_a2)y £A (x-y) (3.35) 

x£A yGZ d \x 

A computation shows that fi(x) — fi(y) = fi(x)[l — go{x — y)] + fo{x)g\{x — y), with go(x) = 
{1,0,-1,0} and gi(x) — go(x — 1). Then we can rewrite (3.35) as 

w i E E e _< ^ L 't /l ^ 1 ^ 1_so ^ xl ~ ! ' 1 ^ + ^ ^ xl ^ £ ' 1 ^ xl ~ yi ^ (— l) x ^y+( ;ri "fi)(_ 9 2 )y eA (x — y) (3.36) 
xeAyez d \x 
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Performing the summation over x yields 



^ 2 |A| J2 cos (f(l - go(xi))) cos fli^i)) (-l) x+ail (-^)n A (x) = 



Pi 



= vf\A\ [co - sin 2 . 9o (zi))(-l) x+a:i (-d 2 )y £A (x) 

x^O 

The conclusion is 

-L 2 ( U i)*W^ 1 (x-y)^ = asin 2 § 



|A| 



(3.37) 



(3.38) 



x,y£A 



where a = X) x ^o(l — .9o(^i))( — l) x+a:i 5 2 F eA (x). We will show in Appendix A that a > 0. Let 
us now consider yeA (M x )*W22(x — y)w y - This double summation is equal to 

^E E e- i *W 1 ^ 1 )-/i(»0)/ ( a;i )/ ( 1 , 1 )(_l)(x3-»s) + -+(x < ,- W ) ( _^ ) y eA(x _ y) (339) 

xeAyez d \x 

Using that fi(x)f 1 (y) = g (x-y) + (-l) x g 1 (x-y) and fo{x)f {y) = g {x-y)- (-l) x g^x - y) , 
we can rewrite this as: 

"a E E e- i *{ /l(a;i)[1 - So(:ci - ! ' l)I+/o(3:i)si(3:i - ! ' l) } [so^i - - (-l) xi gi ( Xl - yi )] ■ 

x£A yGZ d \x 

.(-l)^-y3)+~Mxd-y<i)(_ d ^Y SA (x-y) (3.40) 



Performing summation over x yields 



- 2 2 iaie 

x^O 



g (xi) cos (y (1 - ffo (zi))) cos (yffi(a;i) 
-gi(xi) sin (y (1 - go(xi)j) sin (y#i(a;i)) 



(3.41) 



(-1) 



x 3 H hr d / n2 



(-^)n A (x) 



This can be rewritten as 



« 2 |A| [eo - cos 2 y EC 1 - 0o(zi))(-l) x+X2 (-9 2 )n A (x) 

x^O 



The conclusion is 



jX[ E (-x)^ 2 ° 2 (x-y) U 2 ^ 7 cos 2 ^ 



(3.42) 



(3.43) 



x,y£A 



where 7 = Ex#o( 1 ~ 9o(xi))(-l) x+X2 d$Y £A (x.). We will show in Appendix A that 7 > 0. 
Finally (3.43) is valid even if in the l.h.s we exchange the index 2 with j > 3. Substituting all 
this into (3.32) yields: 



W°(p) > 



/asin 2 ^ 
7 cos 2 ^ 



\ 







V 







7Cos 2 ^-/ 



(3.44) 



By interchanging the roles of p\ , . . . , pa we finally get 

W° (p) > Y (p) , Yij (p) = (a sin 2 | + 7 E cos2 f ) 



(3.45) 
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This bound proves Conjecture 7.9 in [10]. 

Now, the proof of existence of long range order, given the bounds (3.30) and (3.45), is 
standard. We reproduce it here, for completeness. Given I e {1, . . . , d}, we note that, by the 
rotation symmetry of (-)^ A , 

1 1 p 

Therefore, 

]xr^-w < W = 5-]xrEo«(p)- ( 3 - 4? ) 

pjtire 

Let us now consider the sum in the r.h.s. Note that, for any p ^ S, we can use (3.30) and 
(3.45) to conclude that 

Qate) < ± ( q sin2 I + cos2 |) ' • ( 3 - 48 ) 

On the other hand, for p = 7r TO , m / f, we can note that W$(7r m ) = (5^(1 — ^ TO )(ei — e ), 
with ei = Ex#o( _1 ) X+:E2 ( _9 i) y eA( x ) and e o defined in (3.21). Using (3.45), we see that 
ei — eo > (a + -f)/d. Therefore, if v % = Su/(ei — eo), we can use (3.29) to get 

VV [w*(*r m )Q(7r m )l^(7r m )l v* = Qu(n m ) < < (3.49) 

z — ' L J ij Zp ei — eo 2p a + 7 

and we conclude that (3.48) is valid for p = 7T TO , m ^ f , as well. Substituting (3.48) we get 

— (Si t St„ t ) PtA > - - ^ j£ asin 2 f + 7 E^cos2f ' (3 ' 50) 



If we note that 



<£,#-w,W = m E (3-5i) 

1 1 x.yGA 



and take the thermodynamic limit in (3.50) we finally get (2.5), with 



Cdi(3) 1 2/3 (2^ Q sin 2 f + 7 £^ cos 2 f ' (3 " 52) 
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APPENDIX A: THE CONSTANTS a AND 7 ARE POSITIVE 

In this Appendix we show that the constants a and 7 in (3.45), defined right after (3.38) and 
(3.43) respectively, are positive. Let us first consider a = X)x^o(l — 9o( x i))(—]-) x+Xl diY £A (x) . 
Note that 1 — go(xi) is a non-negative even function of x\, vanishing at x\ = 0. Using (3.12) 
we can rewrite 

a = TT^rT E (i-SoOsi)) / + 4 e ^x+*i]-xx e -* lV ^+i*; f (A-1) 
y Z7T ) x 1>0 J 
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where ttI is the {d — l)-dimcnsional vector with components all equal to it. The summa- 
tion over xj_ can be explicitly computed and produces a (d — l)-dimcnsional delta function 
(27:) d - 1 Emez*- 1 S ( k ^- + n ( 2m + !))■ So we find 



a=J2 £ (! - 9o(x 1 )) y Jn 2 {2m. + l) 2 + e 2 e -*iV* 2 (2™+i) 2 +^ > (A.2) 

and the proof that a > is concluded. 

Similarly, let us consider 7 = £ x#0 (l - .9o(a;i))(-l) x+:E2 <9fy £A (x). Note that - 
go(xi)) is an even function of x\, vanishing at x\ = with alternating signs. Using (3.12) we 
can rewrite 

^=7^t E (-ir +1 (i-aa(»i))/ 7^= r ^-e-^V^ 

(A.3) 

Performing the summation over xj_ and using that (— l) Xl+1 (l — goO^i)) = {0, 1, —2, 1}, we find 



ATT 2 m\ cosh yjir 2 (2xn. + ttT - <§i) 2 + e\ - 1 

7 ~ V^ 2 (2m + 7rl-e 1 )2 +e 2 shih2^2^+ ttI - e^ 2 + ef 

where ei = (1,0,..., 0). The proof that 7 > is concluded. 
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